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ABSTRACT 

We report the results of a wide-field narrowband survey for redshift z ~ 5.7 Lya emitters carried 
out with the 34' x 27' field-of-view SuprimeCam mosaic CCD camera on the Subaru 8.3-m telescope. 
Deep narrowband imaging of the SSA22 field through a 120 A bandpass filter centered at a nominal 
wavelength of 8150 A was combined with deep multicolor RIz' broadband imaging with SuprimeCam, 
then supplemented with BVRZ imaging taken with the 42' x 28' field-of-view CFH12K camera on 
the Canada-France-Hawaii 3.6-m telescope to select high-redshift galaxy candidates. Spectroscopic 
observations were made using the new wide-field multi-object DEIMOS spectrograph on the 10-m 
Keck II Telescope for 22 of the 26 candidate objects. Eighteen of these objects were identified as 
z ~ 5.7 Lyman alpha emitters, and a further nineteenth object from the candidate list was identified 
based on an LRIS spectrum. At the 3.3 A resolution of the DEIMOS observations the asymmetric 
profile for Lya emission with its steep blue fall-off can be clearly seen in the spectra of the identified 
galaxies. This is by far the largest spectroscopic sample of galaxies at these redshifts, and we use it 
to describe the distribution of equivalent widths and continuum color break properties for identified 
Lya galaxies compared with the foreground population. The large majority (at least 75%) of the 
lines have rest frame Lya equivalent widths substantially less than 240 A and can be understood in 
terms of young star forming galaxies with a Salpeter initial mass function for the stars. With the 
narrowband selection criteria of (I — N) > 0.7 and N < 25.05 (AB mags) we find a surface density of 
Lya emitters of 0.03 per square arcminute in the filter bandpass (Az = 0.1) down to a limiting flux of 
just under 2 x 10~ 17 erg cm -2 s — 1 . The luminosity function of the Lya emitters is similar to that at 
lower redshifts to the lowest measurable luminosity of 10 43 erg s" 1 as is the universal star formation 
rate based on their continuum properties. However, we note that the objects are highly structured 
in both their spatial and spectral properties on the angular scale of the fields (w 60 Mpc), and that 
multiple fields will have to be averaged to accurately measure their ensemble properties. 

Subject headings: cosmology: observations — early universe — galaxies: distances and redshifts - 
galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The study of high-redshift galaxies in the early universe 
has seen substantial progress in recent times, with the 
discovery of objects out to z = 6.5 (Hu et al. 2002, Ko- 
daira et al. 2003) and with the rapidly incre asing num 
berof galaxies discovered at redshifts z » 5 ijD ev et al. 
19981 Hu, Cowie, & McMahon 1998: iWevmann et al " 
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with radio galaxies (|van Breugel et all 11999ft and AGN 
discovered in deep X-ray exposures with CHANDRA 
l)Barger et al.l l20021 have been seen out to z ~ 5.19, and 
at optical wavelengths quasars have been found with the 
Sloan survey out to z = 6.42 ijFan et al.ll200ll [20031. 
However, because of the relatively sparse distribution of 
detectable z > 5 galaxies it has hitherto been difficult to 
confirm more than a few at a time, and thus to develop 
the large samples needed to understand the robustness of 
candidate selection criteria and the spatial structure of 
objects on the sky. Consequently there are large uncer- 
tainties in the inferred luminosity function and universal 
star formation rates. 

The availability of mosaic CCDs allows wide-field cov- 
erage of the sky on scales of half a degree to a degree 
at a time (e.g., [Malhotra fc Rhoadsl 120021 lOuchi et alJ 
l2003trAiiki et al.l l2003). This gives us the opportunity to 
generate large samples of objects to examine their sta- 
tistical properties and construct more precise luminosity 
functions. 

This paper is the first in a series which will describe 
searches for Lya emission- line galaxies at z = 5.7 and 
z = 6.5 in a number of such wide fields. It repre- 
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sents a continuation of our narrowband filter survey 
llCowiefcHull7993 IHn et alJlTflflSl OH |H|, which 
probed for Lya emission from z — 3.4 — > 6.5 using 
smaller area s urveys. We use the wi de-field capabilities of 
SuprimeCam l)Mivazaki et all2002ft on the Subaru 8.3-m 
Telescope to obtain samples of narro wband selected Lya 
emitters. Following IHu et all l)1998j) we use narrowband 
filters corresponding to z ~ 5.7 and z ~ 6.5 Lya emis- 
sion lines which lie in dark spectral portions of the night 
sky. We combine these with extremely deep continuum 
observations to select the high-z emission-line galaxies 
in a statistically uniform fashion. We then use the new 
wide-field DEIMOS spectrograph on the Keck II 10-m 
telescope to obtain a nearly complete spectroscopic fol- 
lowup. 

In the present paper we describe our survey of z ~ 5.7 
Lyman alpha emitters i n a 700 square a rcminute area 
around the SSA22 field i jCowie et al.lll994|) . The exis- 
tence of extremely deep multi-color BVRIz' continuum 
images for this field allows a detailed examination of the 
effectiveness of color selection criteria in identifying Lya 
candidates a t z ~ 5.7, similar to our studies at z ~ 3.4 
iCowie fc Hi (|1998t) . In particular, we can examine the 
robustness of color thresholds for the selection of Lya 
emitters as contrasted with [OIII] and [Oil] emission- 
line galaxies. 

Discriminating Lya emitters from lower z emission- 
line objects can be difficult at lower redshifts 
l)Stern fc Soinradlll999t) . but at higher redshifts the con- 
tinuum break signature becomes so extreme that there is 
little likelihood of a misidentification provided we have 
sufficiently deep continuum images. Songaila and Cowie 
( 2002) have used spectra of the bright quasars to mea- 
sure the average transmission of the Lya forest region as 
a function of redshift over the z = 4 — > 6 range. Their 
results translate to a magnitude break across the Lya 
emission line of 



Am = 3.8 + 20.3 logi 




For z = 5.7 systems this corresponds to a 3.4 
magnitude break, and color measurements alone can 
be used to provide a substantial selection of high-z 
galax y candidates (e. g ..llwata et a Lll2003t I Stanwav et al.l 
200l IBunker et alJ 120031 ILehnert k Bremer] I2003t 
Bouwens et al .1120031: [Dickinson et alJl2003|) . However, 
because the continuum magnitudes are extremely faint 
only the deepest observations currently feasible can find 
such objects, while the emission-line objects are much 
easier to find. Furthermore, continuum selected objects 
can be extremely hard to confirm spectroscopically, while 
the narrowband selected objects are generally easy to 
identify. Thus the narrowband filter samples provide an 
invaluable tool for identifying large homogeneous sam- 
ples of high-redshift galaxies. 

With high-resolution spectroscopic followups of the 
candidates it is possible to both confirm that the emis- 
sion is due to red-shifted Lya and to estimate the statis- 
tics of interlopers. For SSA22 our narrowband and 
color selections yield an initial sample of 26 candidate 
z ~ 5.7 emitters in the field. The powerful new wide-field 
spectroscopic capa bilities of the DEIMOS spectrograph 
ijFaber et alJ2003T) used on the Keck 10-m telescope, with 
an approximate field of 16.'7 x 5.'0 and good red sensitiv- 



ity, allow us to follow up these candidates in an extremely 
efficient manner. We have now observed 23 of the 26 ob- 
jects with either this instr ument or with the LRIS spec- 
trograph l|Oke et alJll995j) on Keck. This results in the 
confirmation of 19 new redshift z ~ 5.7 galaxies. The 18 
candidates observed with DEIMOS unambiguously show 
asymmetric profiles with the high resolution observing 
mode possible with this red-optimized spectrograph. 

In Section [21 we summarize the narrowband and con- 
tinuum imaging together with the spectroscopic observa- 
tions. Section describes the color selection and color- 
color plots for our selection criteria. The spectroscopic 
identifications and colors of candidates are tabulated and 
labeled in these diagrams. In Section 0] we discuss the 
spectroscopic identifications of Lya emitters, and the in- 
cidence and properties of foreground emitters. In Sec- 
tion [S] we discuss the spatial structure of the objects in 
the fields and show that they are highly correlated. We 
compare the shapes and equivalent widths of the lines 
with those of lower redshift samples and discuss their 
interpretation. Finally, the Lya and UV continuum lu- 
minosity functions of the identified z ~ 5.7 candidates 
are constructed and compared with lower redshift sam- 
ples. We argue that the star formation rates at z ~ 5.7 
are probably at least comparable to those at lower red- 
shifts, but postpone a more detailed discussion until we 
can use an ensemble average of a number of wide fields. 

2. OBSERVATIONS 
2.1. Optical Narrow Band- Selected Survey 

The present images were obtained using a 120 
A (FWHM) filter centered on a nominal wave- 
length of 8150 A in a region of low sky back- 
ground between OH bands. (The nominal spec- 
ifications for the Subaru filters may be found at 
http : //www . naoj . org/Observing/Instruments/SCam/s 
and are also described in (Aiiki et al. 2003).) This cor- 
responds to Lya at a redshift z ~ 5.7, similar to the 
narrowband wavelengths probed in earlier imaging 
studies wit h the LRIS ins trument on the Keck 10-m 
telescopes l|Hu et al.l |l998). The filter profile in the 
//1.86 beam is shown in Figure ^ The Gaussian 
shape of the filter profile used on SuprimeCam may 
be compared with the squarer profile of the 108 A 
bandpass 8185 A filt er used in the L RIS parallel beam 
for the Keck studies l|Hu et alJ ll999). The narrowband 
NB816 filter is well-centered in the Cousins I band 
filter, which we use as a reference continuum bandpass 
for our detection of Lyman alpha emission. Hereafter, 
we refer interchangeably to this filter as NB816 or as N. 
All narrowband and Z-band magnitudes are given in 
the AB system. 

Narrow-band images were taken with SuprimeCam on 
the 8.3-m Subaru telescope on the nights of UT 2001 
25-26 June and 2001 20 October under photometric or 
near photometric conditions. The imaging observations 
are summarized in Tabled with exposures made during 
times of heavy cirrus shown in parentheses. The data 
were taken as a sequence of dithered background-limited 
exposures and successive mosaic sequences were rotated 
by 90 degrees. Corresponding continuum exposures were 
always obtained in the same observing run as the narrow- 
band exposures to avoid false identifications of transients 



LUMINOSITY FUNCTION OF REDSHIFT z ~ 5.7 LYMAN ALPHA EMITTERS 3 

Table 2. Summary Log of CFH12K SSA22 Imaging 




"10 ; i i i i I i i i i I i i i i I i i i i I i i i i I i i i i : 

6500 7000 7500 8000 8500 9000 9500 

OBSERVED WAVELENGTH (A) 

Fig. 1. — Plot of the NB816 and Cousins J filter profiles in the 
f/1.86 beam of SuprimeCam superposed on a night-sky spectrum 
obtained with the LRIS spectrograph on Keck. The narrowband 
filter is situated in a region of low sky background and is well 
centered in the /-band filter. Due to the f/1.86 beam, the profile 
has a Gaussian shape, in contrast to the squarer profiles that would 
apply for parallel or slower beams. Peak throughput is ~ 85%. The 
strength of night-sky OH airglow lines is highly variable with time, 
and the spectrum is meant to illustrate regions of high and low sky 
background rather than providing an absolute calibration. 



Table 1. Summary Log of SuprimeCam SSA22 Imaging 
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such as high-z supernovae, or Kuiper belt objects, as 
Lya candidates. A detailed descript ion of the full reduc - 
tion procedure for images is given in lCapak et alJ |2004). 
In addition to the SuprimeCam multi-color imaging on 
SSA 22, wide-field observat ions with the CFH12K cam- 
era ijCuillandre et al]l200(1) were obtained (Table EJ for 
a total of 5hr 20min in B, 4hr 15min in V, 7hr 3min in 
R, and 7hr 15min in Z on the CFHT 3.6-m telescope. 

The SuprimeCam data were calibrated using the 
photometric and spectrophoto metric standard, HZ4 
ijTurnshek et alll990UOkell990l). and faint La ndolt stan- 
dard stars in the SA 95-42 field <jLandoltll99'2T) . An astro- 
metric solution was obtained using stars from the USNO 
survey. Limiting magnitudes (5(7 for 3" diameter aper- 
tures, expressed as AB mags) are: 25.5 (B), 25.6 (V), 
26.1 (R), 25.7 (I), 25.0 (z'), and 25.2 (NB816). 

2.2. Spectroscopy 

We used the new wide-field DEIMOS spectrograph on 
the Keck 2 10-m telescope for spectroscopic followups of 
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candidate z ~ 5.7 Lya emitters. The G830 £/mm grat- 
ing blazed at 8640 A was used with 1" wide slitlcts. In 
this configuration, the resolution is 3.3 A, which is suffi- 
cient to distinguish the z ~ 1.19 [Oil] doublet structure 
from the profile of redshifted Lya emission (Fig.|2J. Red- 
shifted [OIII] emitters (z ~ 0.62), show up frequently as 
emission-line objects in the narrowband and can easily be 
identified by the doublet signature. The observed wave- 
length range (AA ~ 3840 A, with typical coverage from 
~ 5900 — 9700 A) generally also encompasses redshifted 
H/3 and [O III] lines for instances where emission-line can- 
didates might be Ha at z ~ 0.24, and deals with the 
problematic instance of extragalactic H II regions with 
strongly suppressed [Nil]. The G830 grating used with 
the OG550 blocker gives a throughput greater than 20% 
for most of this range, and - 28% at 8150 A. The obser- 
vations were taken on UT 12-13 September 2002, with 
the first night showing light to moderate cirrus, and the 
second night clear, and on two essentially clear nights on 
UT 28-29 August 2003. 

Earlier spectra of a preliminary cut of candidates were 
also obtained with the LRIS spectrograph during runs in 
UT 18-19 August 2001 and 21 October 2001. The ob- 
servations were made in multislit mode using 1.4" wide 
slits with the 400 f/mm grating blazed at 8500 A, which 
in combination with the slit width gives a resolution of 
roughly 11.8 A. Because of the much higher red through- 
put of the DEIMOS CCDs and optics, and its wider field 
coverage, nearly all of the spectroscopic identifications 
use the DEIMOS reduced data. 

3. CANDIDATE SELECTION 

Figure 13 shows the / — Nab color as a function of nar- 
rowband (Nab) magnitude. All magnitudes were mea- 
sured in 3" diameter apertures, and had average aper- 
ture corrections applied to give total magnitudes. To a 
narrowband magnitude of 25.05, conservatively brighter 
than a 5 a cut of 25.2, there are 127 objects with I — Nab 
colors > 0.7 over the central 26.'5 x 26.'5 field, which was 
conservatively chosen as a region of uniform mosaic cov- 
erage, free of any possible edge-of-ficld illumination vari- 
ation effects. This is about 0.4% of the 38451 objects 
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Fig. 2. — Composite DEIMOS spectra taken through 1" wide 
slitlets with the G830 grating of [OIII] (top panel), [Oil] (middle 
panel), and Lya (bottom panel) emission-line objects lying within 
the NB816 bandpass. At the 3.3 A resolution of this configura- 
tion these systems are easily distinguished. Dashed vertical lines 
mark the positions of the two components of the oxygen doublet. 
The stacked composite profile of the Lya emission systems clearly 
show the steep drop-off on the blue side of the emission feature. 
The superposed dashed profile shows the instrumental response to 
neighboring night sky lines for comparison, and indicates the blue 
side of the emission feature is consistent with an abrupt cut-off. 
The asymmetry of the profile for the Lya emitters is a clear signa- 
ture of the forest absorption to the blue of the Lya emission line. 



with Nab < 25.05 in this area. As can be seen from 
Fig. |21 a simple emission criterion chooses out a mixture 
of [Oil] and [OIII] emitters in addition to the z = 5.7 
objects and a small number of stars. The emission-line 
excess objects have been labeled on this plot according 
to their subsequent spectroscopic identifications: filled 
squares for Lya emitters, triangles for [O II] emitters, di- 
amonds for [O III] emitters. 

The high-z objects are also expected to have very red 
colors because of the effects of the Lya forest scatter- 
ing and of the intrinsic Lyman continuum breaks in the 
galaxies themselves. For a z ~ 5.7 object whose Lya 
emission is in the narrow band, the scattering in the R 
band will be a mixture of Ly/3 from the intergalactic 
medium at z ~ 5.4 and Lya from z ~ 4.5. Their com- 
bined effects result in a transmission of only about 20% 
of the light at this wavelength based on measurements 
of brighter quasars (Songaila 2004, in preparation). This 
corresponds to an R — z' AB break of about 1.8 mags so 
that we expect the actual R — z' colors for the objects 
to be redder than this value. By contrast we expect that 
[Oil] lines which have the Balmer and 4000 A breaks 
between the two bands will have weaker R — z' AB breaks 
of up to about a magnitude, while [O III] or Ha lines 
should have nearly flat continua over the R — z' wave- 
length range. At wavelengths less than about 6000 A we 
expect the galaxies to be faint owing to Lyman contin- 
uum absorption and that this should result in a z ~ 5.7 
galaxy being undetected in the B and V bands. 

We illustrate this in Figure0]where we show the R— z' 
color versus the narrow band excess I — N for all the 
brighter objects with Nab < 23.8 in the field. Known 
[O II] and [O III] emitters are marked. At this bright nar- 
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Fig. 3. — Plot of emission-line excess objects (I — N vs. N) found 
in the NB816 narrowband filter as a function of object magnitude. 
Narrow-band m agnitudes are given in the AB system llOke fc Gunnl 
1983; Okc 1990) of equivalent visual magnitudes. The distribution 
of field objects is also shown as a scatter plot. Filled symbols show 
the properties of spectroscopically identified emission-line objects 
(triangles are [Oil] emitters, diamonds are [OIII] emitters, and 
squares are candidate z ~ 5.7 Lya emitters) and crosses inside open 
diamonds are candidates which were spectroscopically observed but 
not identified. An (I — N) > 0.7 and N < 25.05 selection criterion 
(dashed lines) picks out 127 objects which are shown with crosses 
if they have not been spectroscopically observed. Asterisks mark 
the positions of known stars. 



rowband magnitude cut the color regions populated by 
stars and foreground galaxies are very well defined, and 
we can use this to define a region (solid lines) in which 
z ~ 5.7 Lyman alpha galaxies can be unambiguously dis- 
tinguished from foreground populations. The regions of 
the [O III] and [O II] emitters can clearly be seen. As ex- 
pected, the [O III] emitters are nearly flat in R—z' AB with 
typical values of 0.0 to 0.3. These objects can have ex- 
tremely high equivalent widths, however. By contrast the 
[O II] emitters have lower equivalent widths, with only a 
few exceeding I — Nab ~ 0.7, but have redder colors. 
The higher equivalent width [O II] objects in Figure 0] 
appear to be slighter bluer than the weaker ones, with 
typical R — z' AB colors of just under one. At these mag- 
nitudes there are no objects in the field with red enough 
colors to be a z ~ 5.7 emitter with / — Nab > 0.7. 

Figure 0] shows that red stars can lie in a simple 
I — Nab > 0.7 plus red color break selected sample 
since at the red end the star track angles into this re- 
gion. The reason for this is that the narrow band lies 
between strong stellar absorption features lying in the / 
band and these strengthen for redder stars through to L 
types. This results in an apparent excess in the narrow 
band. This result can be used beneficially in color selec- 
tion procedures to separate out stars from high-z galaxies 
or quasars (Kakazu et al. 2004, in preparation) but for 
the present work we must also eliminate the star track. 
We have therefore used a slightly higher / — Nab bound- 
ary at the reddest R — z' AB levels which is illustrated in 
Figure 0] 

The stars can also be distinguished from candidate Ly 
a emitters on the basis of their Nab — z' colors. The lo- 
cation of the star track is shown by the I — z'vs Nab — z 1 
color-color plot of Fig. [3] The dashed line gives the 
I — Nab > 0.7 selection criterion, and the stars and 
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Fig. 4. — Plot of emission-line excess objects (R — z 1 vs. / — AT) 
found in the NB816 narrowband filter for a bright N < 23.8 se- 
lection. The distribution of field objects is also shown as a scatter 
plot. Filled symbols show the properties of spectroscopically iden- 
tified emission-line objects (triangles are [Oil] emitters, diamonds 
are [O III] emitters, and open diamonds are candidates which were 
spectroscopically observed but not identified. Plus signs indicate 
objects which satisfy the (I — N) > 0.7 condition (shown as the 
dashed line), but which have not been spectroscopically observed. 
Asterisks show spectroscopically identified stars. The adopted se- 
lection criteria for unambiguous z ~ 5.7 Lya emitters in this color- 
color plane is shown by the solid lines. This figure may be compared 
to a similar color plot (Fig. [HJ for a faint narrowband selection, 
23.8 < N < 25.05, which spans the magnitude range of the newly 
discovered z ~ 5.7 Lya galaxies. 



boxes show the track for low-mass stars (including L and 
T dwarfs) for this filter system. Because the Lya emis- 
sion falls within the / band and can be a substantial con- 
tributor to the measured broad-band flux, these objects 
can have blue I — z' colors. In general for / — Nab > 0.7 
selected objects the emitters lie at Nab — z' < 0.2 and 
the stars lie at Nab — z' > 0.2. Because the level of star 
contamination is low and we wished to avoid eliminating 
ambiguous cases we did not apply this criterion prior to 
the spectroscopic observations and chose instead to leave 
a small number of stars in the candidate sample. We will 
discuss these objects in the section on the spectroscopy. 

Our final sample of objects was then taken as all ob- 
jects with Nab < 25.05 which clearly lay within the des- 
ignated area of R — z', I — Nab color-color space and 
which were also not detected in V or B, together with ob- 
jects which satisfied the emission line criterion and were 
not detected at any wavelength shortward of I. Their 
R — z' vs. / — Nab color properties are shown in Fig. El 
with object symbols as previously defined. Arrows show 
2 a limits on the colors. 

This final candidate sample contains 26 objects. Tabu- 
lated coordinates, multi-color magnitudes, and redshifts 
(where measured) for these objects are summarized in 
Table [3J Thumbnail finding charts in the narrowband 
(Fig. EJ) and i?-band (Fig. |HJ are shown, with the ob- 
jects encircled. Each thumbnail is 30" on a side, with 
the circles 2" in radius. 

4. SPECTROSCOPIC CLASSIFICATION 

Twenty-two of the 26 objects in the candidate list were 
observed with DEIMOS with exposure times ranging 
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Fig. 5. — I — z' vs. N — z' color distribution of spectroscopically 
identified narrowband excess objects with I — N > 0.7 (dashed 
line) compared with the spectroscopically identified stars in the 
field. Object identifications have been indicated with the usual plot 
symbols. Squares are Lya emitters, diamonds are [O III] emitters, 
triangles are [O II] emitters and asterisks are stars. The distri- 
bution of L and T dwarf stars is shown by the small open boxes 
which are drawn from the literature (Kakazu et al. 2004, in prepa- 
ration) and lie at redder N — z' colors (dotted line) compared to 
the identified Lya emitters. Two objects in TableEJwhich are spec- 
troscopically unidentified but which are probably stars are marked 
with open diamonds. 




I-N (AB) 

Fig. 6. — Color distribution of narrowband excess objects in 
R — z' for the N > 23.8 and N < 25.05 sample. These data may 
be compared with the bright (N < 23.8) sample color distribution 
of Fig. [I] where the color tracks of galaxies and stars and the loca- 
tions of foreground emission-line galaxies may be more clearly seen. 
Object identifications have been indicated with the plot symbols. 
Squares are Lya emitters, triangles are [O II] emitters, diamonds 
are [OIII] emitters, plus signs are (/ — N) > 0.7 (dashed line) ob- 
jects which were not spectroscopically observed, and asterisks are 
stars. The adopted selection criteria for unambiguous z ~ 5.7 Lya 
emitters indicated in Fig. |l]is shown by the solid lines. R — z' 
color breaks for faint objects are generally limited by the R band 
non-detection sensitivity. 



from 1.5 to 4.5 hours. A randomly chosen sample of 21 
additional objects with (/ — N) > 0.7 which did not sat- 
isfy the color selection was also observed with the choice 
of objects in this case being determined by the mask con- 
figurations. One additional object in the candidate list 
which had been previously observed with LRIS was not 
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Table 3. Properties of Candidate Sample 



Object 


RA(2000) 


Dcc(2000) 


N(AB) 


Z(AB) 


J 


R 


V 


B 


z 


I 


334 220001 


277625 


25.0 


28.7 


28.6 


99.0 


99.0 


99.0 


5.659 


■ 


334 22Q004 


093959 


23.9 


26.5 


26.0 


99.0 


99.0 


27.7 


5.676 


3 


334 234985 


246227 


24.5 


24.9 


26.3 


28.3 


99.0 


99.0 


5.660 


4 


334 27301 


21 6856 


24.4 


25.8 


26.1 


99.0 


99.0 


29.4 


5.670 


5 


334 278015 


0.206272 


24.7 


25.2 


26.3 


28.3 


99.0 


99.0 


5.651 


5 


334 279999 


462454 

. TtUiiTttjTt 


24.6 


25.6 


26.6 


30.4 


99.0 


27.9 




7 


334 201 002 


21 1 263 


24.6 


24.5 


25.6 


99.0 


99.0 


31.4 


5.626 


g 


334 308014 


316606 


24.9 


24.5 


25.6 


26.9 


99.0 


99.0 




9 


334 31 70Q3 


223775 


24.6 


25.2 


25.8 


28.3 


99.0 


29.3 


5.645 


10 


334 320007 


359441 


25.0 


25.5 


25.8 


99.0 


99.0 


99.0 




11 


334 337006 


335259 


24.5 


27.0 


26.2 


27.8 


27.0 


27.2 


5.669 


12 


334 337006 


293762 


24.1 


24.8 


25.6 


99.0 


99.0 


99.0 


5.667 


13 


334 351000 


197130 


24.7 


27.3 


27.9 


28.4 


26.9 


99.0 




14 


334 360005 


321 620 


24.2 


25.0 


25.2 


27.1 


28.3 


28.2 


5.641 


15 


334 381 080 


1 601 89 


24.5 


25.9 


26.5 


28.6 


99.0 


27.9 


5.684 


16 


334.388000 


0.371196 


24.3 


26.1 


25.9 


28.0 


26.6 


99.0 


5.653 


17 


334.389008 


0.390323 


24.6 


23.9 


25.5 


27.6 


99.0 


26.8 


star? 


18 


334.415985 


0.307220 


24.8 


26.6 


26.4 


28.4 


28.1 


26.7 


5.723 


19 


334.420013 


0.404125 


24.9 


25.7 


25.7 


28.8 


27.8 


99.0 


5.633 


20 


334.420013 


0.118385 


24.5 


25.4 


26.6 


99.0 


29.2 


99.0 


1.173 


21 


334.425995 


0.128724 


24.4 


24.9 


26.1 


99.0 


28.3 


99.0 


5.687 a 


22 


334.492004 


0.228247 


24.9 


25.0 


25.7 


27.9 


99.0 


99.0 


0.492? 


23 


334.509003 


0.242099 


24.1 


25.9 


25.9 


99.0 


99.0 


99.0 


5.673 


21 


334.519012 


0.239439 


24.8 


27.4 


28.7 


99.0 


99.0 


99.0 


5.680 


25 


334.548004 


0.083544 


23.8 


25.1 


25.9 


28.1 


99.0 


28.4 


5.709 


26 


334.556000 


0.095670 


25.0 


27.2 


26.9 


99.0 


99.0 


99.0 


5.682 



Note. — Magnitudes are measured in 3" diameter apertures. An entry of '99' indicates that no excess flux was measured. 
i Based on lower resolution LRIS data. 
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Fig. 7. — Narrowband 8150/120 A images of the 26 emission-line candidates selected in Table 151 Each field is 30" on a side. A 2" 
radius circle identifies each candidate, which is labeled by object number and redshift (for confirmed high-2 galaxies). This figure may be 
contrasted with the corresponding i?-band images of Fig. 181 



included in the DEIMOS masks. The DEIMOS masks 
were filled with color selected and magnitude selected 
samples which will be described elsewhere and contained 
just over 900 objects. 

All of the spectra (emission line and field objects) were 
spectroscopically classified without reference to their nar- 
row band strengths and color properties. As can be 



seen from Fig. at the DEIMOS resolution [Oil] and 
[OIII] lines are easily distinguished based on their dou- 
blet structure, while Ha emitters can usually be picked 
out on the basis of the [N II] or [S II] doublets or using 
the H/3, [OIII] complex. Only a small number of emis- 
sion line objects show single broad asymmetric lines of 
the type illustrated in Fig. |3 and these were classified as 
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Fig. 8. — R-band images of the selected emission-line candidates shown in Fig.|7| A 2" radius circle indicates the candidate's location in 
each 30" X 30" field. Candidates are labeled by object number and redshift (for confirmed high-z galaxies) from Table 151 Absorption due 
to the Lyman alpha forest should suppress ij-band light from redshift z ~ 5.7 galaxies. 



Lya emitters. Where there was uncertainty in the classi- 
fication the object was marked as observed but unidenti- 
fied. Where an object was observed multiple times each 
spectrum was measured individually. However, there was 
no disagreement among identifications and in the final 
compilation these spectra were averaged. 

The redshift identifications of the high z candidate list 
are summarized in Tabled! Of the 22 DEIMOS observed 
objects 18 were identified as Lyman alpha emitters, as 
was the LRIS observed object. One object was classi- 
fied as an [O II] emitter and three could not be clearly 
identified. The spectra of the identified (N — I) > 0.7 
objects are shown in Figures |H1 El an d EJ O ne of the 
unidentified objects has a weak continuum spectra which 
probably corresponds to a star and lies in the area of 
z 1 — N, I — N color-color space consistent with this in- 
terpretation as does one of the unobserved objects (17 
and 8 respectively). One object (21) has a strong red 
emission line which is symmetric and narrow and which, 
if it is Ha, places the object at redshift z — 0.492. There 
is a weak feature at the [OIII] A5007 emission position 
which could be consistent with this interpretation but, 
irrespective, this object is clearly not a z ~ 5.7 Lyman 
alpha emitter. Finally, a retrospective analysis of the 
[O II] emission line classified object (20) suggests it is 
ambiguous whether it should have been classified in this 
way since the [O II] lines are slightly displaced from the 
correct positions. It is possible this is actually a Lyman 
alpha emitter. As we show in Fig. 1121 if this object is an 
[O II] emitter its colors are extremely anomalous. 

For the remainder of the work we will focus on the 
19 spectroscopic Lya emitters recognizing that a further 
4 objects could fall in the sample. This small level of 
incompleteness does not affect the discussion in any sta- 
tistically significant way. 



5. RESULTS AND DISCUSSION 
5.1. Spatial and spectral distributions 

In Fieure lT^l we show the distribution of the 19 Lyman 
alpha redshifts compared with the nominal transmission 
profile of the filter shown by the dotted line. All but one 
of the 19 objects lie below the redshift corresponding to 
the 81501 central wavelength of the filter. The [OIII] 
emitters also lie in the lower half of the filter and taken 
together these results could suggest that the true on- 
telescope center of the filter is shorter than the measured 
value. However, both could be the result of structuring at 
the respective redshifts. Even if we shift the filter center 
25 A blueward to match the [O III] data (dashed line) 16 
of the 19 Lya emitters still lie in the lower half of the 
filter. The probability of seeing such a biased result is 
only 0.0003 while with the laboratory traced filter profile 
the probability would drop to 2 x 10~ 6 . It appears that 
the sample is extremely structured in velocity space. 

The objects lying in the lower redshift interval also 
appear to be clustered in space, with many of them lying 
in a long elongated filament running diagonally across 
the field (Figure IT4"}) . All of the objects between z = 
5.632 — 5.672 lie in this region. 

We postpone a more detailed discussion of the corre- 
lation function properties until the analysis of the multi 
field sample but the present results emphasize that single 
fields may vary considerably even with areas as large as 
the present ones and that many independent fields are 
required to average out these systematic effects. 

5.2. Line Shapes 

As can be seen from Figure the Lya shapes are re- 
markably uniform, consisting of a slightly peaked profile 
and a trailing red wing. In Figure lT31 we show the stacked 
profiles of the 18 Lya emitters with DEIMOS spectra 
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Fig. 11. — Spectra of emission line objects with the [O III] doublet 
signature. As in Fig. [5] the bottom spectrum shows the positions 
of the strong nightsky lines which may give rise to subtraction 
residuals. 
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8000 8050 8100 8150 8200 8250 8300 
Fig. 9. — Extracted DEIMOS spectra of single-line emission line 
objects with asymmetric profiles in the candidate sample. Spectra 
are labeled with the assigned object number and redshift given in 
Table 3. The bottommost plots show the positions of strong night- 
sky lines for comparison with subtraction residuals (e.g., objects 
18, 23, 24, 25, 12, 4, and 2). 
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Fig. 10. — Spectra of emission line objects with the [Oil] dou- 
blet signature. Candidate object 20, from Table l3l is identified as 
an [Oil] emitter at redshift z = 1.172. The remaining 8 spectra 
show galaxies selected according to [O II] emission excess within 
the narrowband filter. As for Fig. |^] the bottom spectrum shows 
the positions of the strong nightsky lines which may give rise to 
subtraction residuals. 



Fig. 12. — Color-color plots in R — z' vs. I — AT of [Oil] emit- 
ters (triangles) and [OIII] emitters (diamonds). [Oil] emission-line 
galaxies display a shallow continuum break in (R — z') of about 1 
magnitude. The range of emission-line excess (J — N) reflects the 
cut-off of selection criteria at the low excess end, but is generally 
modest. The [OIII] emitters in the narrow-band filter can dis- 
play stronger observed-frame excess emission, reflected in the wider 
range of emission excess in (I — N). The [OIII] galaxies have gen- 
erally flat spectra (R — z') ~ 0, but there is a bluing trend seen 
in the (R — z') colors for galaxies with very strong [OIII] emission 
lines. The one object classified as an [O II] emitter (object #20) in 
the candidate sample lies in the upper right corner of the left hand 
panel. 



compared with the resolution element of the instrument, 
which clearly illustrates this point, but individual profiles 
are also generally similar in shape and width. This type 
of profile is a simple outcome of the blue truncation of 
the spectrum as is illustrated in Figure lTBl where we show 
a simple model where the intrinsic line is modeled as a 
Gaussian and then is truncated below the object redshift 
by the Lyman alpha scattering of the IGM. When this 
profile is convolved through the instrument resolution it 
produces an observed spectrum almost identical to that 
observed as is shown in the lower panel of the figure. 
The required width of the intrinsic line which matches 
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Fig. 13. — The redshift distribution of Lya emission galaxies 
compared to the transmission profile of the narrow band filter (dot- 
ted line). Nearly all of the objects lie at the low redshift end of the 
selection, This result remains true even if we shift the filter center 
25 A to the blue (dashed line). 
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Fig. 14. — The spatial distribution of the 19 Lya emission galax- 
ies is shown by the large solid (z < 5.687) and open (z > 5.687) 
boxes. The remaining objects in the candidate list are shown as 
crosses for unidentified objects, asterisks for possible stars, and the 
small box for the possibly misclassified [Oil] emitter which, if it 
were a Lyman alpha emitter, would lie in the lower redshift inter- 
val. Many of the low z objects lie in the large diagonal filament 
which runs across the field. Most of the objects in this region lie 
in the more restricted redshift interval z = 5.632 — 5.672. Objects 
lying in this redshift region are enclosed in larger open boxes. 
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Fig. 15. — The stacked Lya emission line profile of the eigh- 
teen new z ~ 5.7 galaxies reported in this paper. Individual ex- 
tracted spectra were normalized to their maximum value, converted 
to their rest-frame wavelength scale using the redshift correspond- 
ing to their peak flux, and co-added. The instrumental profile as 
determined from neighboring resolved nightsky lines is overplotted 
on the peak for comparison. The asymmetric profile can be clearly 
seen. 
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Fig. 16. — A simple model profile for the Lya emission feature 
in Fig. 1151 The top profile shows a Gaussian emission-line profile 
with Doppler width a = 200 km s — 1 with the blue half of the line 
(dashed curve) absorbed by the foreground Lya forest. When this 
is convolved with the instrumental profile of the spectrograph, the 
resulting profile (dashed curve) shows a good fit to the observed 
stacked Lya profile (solid curve) in the bottom half of the figure. 



the observed stacked spectrum corresponds to a a of 200 
km s _1 and this value also produces a reasonable fit to 
nearly all the individual line profiles. This is much wider 
than would be expected from the velocity dispersions of 
the parent galaxies at this redshift and almost certainly 
represents velocity broadening of the Lya during the es- 
cape from the galaxy. Given the complexity of this pro- 
cess there is probably relatively little useful information 
in the widths. The truncation also results in the mea- 
sured redshifts, which were set to the peak of the line, 
being about a 100 km s~ x redward of the true redshift, 

5.3. Equivalent Widths 



Since the narrow band lies almost precisely in the cen- 
ter of the / band the truncation effect reduces the / band 
flux and the narrow band flux in very similar amounts 
and we shall assume that the equivalent widths (EWs) of 
the lines do not need to be corrected for this effect. (This 
approximation does depend on the IGM scattering not 
extending substantially redward of the line center which 
could be incorrect.) Because the filter is not rectangu- 
lar in shape we do need to allow for the filter profile in 
computing the EWs. 

The simplest approach is to compute the equivalent 
widths based on the I and N magnitudes. Defining the 
quantity 

r = icr - 4 *^- 7 ) 
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the observed frame equivalent width becomes 



EW 



R-l 

* R 
~ AA 



where <f> is the narrow band filter response normalized 
such that the integral over wavelength is unity and AA is 
the effective width of the / filter, which is approximately 
1400 A. The narrow band filter is often assumed to be 
rectangular in which case <p becom es 1/AA where AA is 
the w idth of the narrow band fe-g.- lRhoads fc Malhotral 
2002) but as can be seen from FigureEDthis is not a very 
good approximation in the present case. 

For very high equivalent width objects such as the Ly- 
man alpha emitters the denominator in this equation be- 
comes uncertain unless the broad band data is very deep 
and this can result in a very large scatter in the measured 
equivalent widths. Therefore in order to check the result 
we also computed the equivalent widths by extrapolating 
from the z' observations assuming that these provide a 
better measure of the true continuum flux. We assumed 
in this case that the spectral energy distribution was flat 
in frequency space giving the formula 



EW 



2 x Norm x [R 2 - 1) 



where 



R 2 = 10 



-0A*(z'-N) 



and Norm is the square of the wavelength ratio of the 
center of the z' filter over that of the narrow band. In 
order to avoid uncertainties from the exact wavelength 
center of the filter we restrict ourselves to wavelengths 
where the nominal filter response is more than 0.67 of 
the peak. 

The rest frame equivalent widths calculated with both 
methods generally agree well with each other though 
there are significant differences in a small number of indi- 
vidual cases. The distribution based on the extrapolated 
fluxes is shown in Fig. IT7| In contrast to the results of 
Malh otra fc Rhoadsl 1)20021 most of the rest frame equiv- 
alent widths have values less than 240 A, which can be 
understood in terms of galaxies with Salpeter initial mass 
functions. (See Malhotra and Rhoads for an extensive 
discussion of this issue.) As discussed in the previous 
paragraph, the equivalent width determination is very 
sensitive to the depth of the continuum band exposures, 
a nd the difference between th e present results and those 
of lMalnotra fc Rhoadsl (2002) may arise from the deeper 
continuum exposures used here. Only four of the objects 
have poorly defined rest frame equivalent widths and are 
shown at nominal values of 240 A in the figure but even 
these objects could well lie at lower values based on the 
uncertainties in the determination. The overall distribu- 
tion of Fig. El is quite similar to the distri butions mea- 
sured at lower redshifts (e.g. Figure 21 of Fuiit a et alJ 
(2003J which summarizes these measurements.) 

5.4. Lyman alpha Luminosity Function 

Because of the high observed frame equivalent widths 
the Lya fluxes are insensitive to the continuum determi- 
nation. However, they do depend on the filter response at 
the emission line wavelength so we have again restricted 
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Fig. 17. — Rest frame equivalent width distribution for the 15 
Lyman alpha emitters lying in the redshift range where the filter 
profile lies above 0.67 times its maximum value. 
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Fig. 18. — The luminosity function based on the Lya emis- 
sion line for the z ~ 5.7 emitters (filled squares) compared to the 
emission line luminosity function at z ~ 3.4 (open squares). The 
vertical dashed line indicates the limits of the sample, which does 
not yet reach the flat part of the z ~ 5.7 luminosity function. 



ourselves to redshifts where the nominal filter response 
is greater than 67% of the peak value. 

Because the volume is simply defined by the selected 
redshift range the luminosity function may be obtained 
by dividing the number of objects in each luminosity bin 
by the volume. The calculated Lya luminosity function 
is shown in Figure 1181 where we compare it with the 
z = 3 luminosity function computed from the results 
of IHu et alJ l|1998D . The observations of IKudritzki et all 
1120001) give a similar value at this redshift. The present 
function samples generally higher luminosities than the 
z = 3 samples but is broadly consistent in number den- 
sity at the lowest l uminosities. It lies slightly above the 
iFuiita et al.l |2003) values at similar luminosities, but 
this probably reflects the higher equivalent width selec- 
tion that was used in that work. 

5.5. Continuum Luminosity Function 
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Fig. 19. — The luminosity function of the UV continuum of z ~ 
5.7 emitters compared with the Steidel et al. 1999 UV luminosity 
functions at redshifts 3 (diamonds) and 4 (triangles). Measured 
points from the Lya emitters are shown as solid boxes; the open 
boxes show assumed values if emitters pick out 20% of the LBG 
sample as Steidel et al. (2000) find to be the case at z ~ 3. 



We may also view the Lyman alpha selection as pick- 
ing out that subset of the Lyman break galaxies (LBGs) 
at this redshifts which have strong line emission. While 
some of the objects are very faint in the z' band, which 
corresponds to a rest frame of about 1400 A, most (16 
of the 19) have well measured z' magnitudes lying in the 
range from 24 to 27. These magnitudes correspond to ab- 
solute continuum magnitudes at 1450 A which are very 
similar to the typical LBG at z = 3 — 4. The use of con- 
tinuum magnitudes to estimate the star formation rates 
avoids the extremely complex problems of the Lyman 
alpha escape process and the uncertainties in the correc- 
tion of the Lyman alpha fluxes for intergalactic scatter- 
ing which are present in the determination of the Lya 
luminosity function. 

Figurc lT^l shows the continuum luminosity function de- 
rived from the distribution of the z' magnitudes. Once 
again we have restricted ourselves to objects lying in the 
higher transmission region in the filter. The solid boxes 
show the measured values with 1 a Poisson uncertain- 
ties based on the number of objects in the bin. There 
are no objects with absolute magnitudes brighter than 
—22 and for the final data point we show a 1 <r upper 
limit. No correction has been applied for dust. Fig- 
ure El also shows the raw (dust uncorrected) UV lu- 
minosity function for Lyman break selected galaxies at 
redshifts z ~ 3 (open diam onds) and z ~ 4 (open tri- 
angles! (iSteidel et all I"l999. 2000). At these redshifts 
ISteidel et alJ l)2000|) found the Lyman alpha selection 
picks out about 20% of the Lyman break galaxies. The 
open squares show the values which would be obtained 
for the LBG luminosity function at z = 5.7 if the Lyman 
alpha selected fraction were similar to that at the lower 
redshift. It matches extremely closely to the lower red- 
shift luminosity functions. The Lya selected fraction at 
these redshifts could be larger, and compariso n with the 
color-selected samples of lStanwav et alJ l)2003|) would be 
consistent with this interpretation. 

Once again, recognizing the systemic uncertainties in a 



single field observation and the uncertainties involved in 
the Lya selected fraction, we postpone a more detailed 
discussion to a subsequent paper. However, it is clear 
already from Figure IT§1 that star formation rates at z = 

5.7 are roughly similar to those at lower redshift and 
there is no rapid decline in the star formation rates at 
these redshifts. 

6. CONCLUSIONS 

We may summarize the major conclusions of the paper 
as follows. 

1) Combined emission line and color break techniques 
provide an extremely efficient method of selecting objects 
at z = 5.7. Selection criteria of (N - 1) > 0.7, (R — z') > 

1.8 and (N - z') < 0.2 pick out 24 objects in the field. 
At least 19 of these are spectroscopically confirmed as 
Lyman alpha emitters at this redshift. Only one or at 
most two of the sample selected in this way are lower 
redshift interlopers. The accuracy is sufficiently good to 
allow use of color selected samples without the necessity 
of spectroscopic followup for many purposes. 

2) Even on scales as large as the present 26' field the 
distribution of objects is highly structured in both the 
spectral and spatial dimensions. Individual fields could 
easily show factors of several variation in the number 
density of objects selected in a 100 A bandpass and a 
substantial number of fields need to be observed to aver- 
age out these effects. 

3) The properties of the Lyman alpha selected galaxies 
are quite similar to those observed at z w 3. 

4) Star formation rates at z = 5.7 are roughly simi- 
lar to those at z = 3 and z = 4 though there are very 
substantial uncertainties in this estimate. 
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